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I. INTRODUCTION 


Two recent NASA earth-orbiting satellites, GEOS-3 and SEASAT-1, 
have carried pulsewidth-1 imited radar altimeters with provision for 
sampling a number of points in the individual radar return waveforms^. 

For these and similar systems, the mean return waveform is [1] the 
convolution of: A) the average impulse response of the quasi-calm sea 
surface, B) the sea surface elevation distribution, and C) the radar 
system point-target response (transmitted pulse as affected by transmitter 
and receiver bandwidths;. The first term. A, includes effects of antenna 
beamwidth and off-nadir pointing angle; "quasi-calm" emp!)asizes that an 
incoherent surface scattering process is assumed but that the sea surface 
elevation distribution is separately written in B in the above convolution. 

A number of papers have described the extraction of ocean significant 
waveheight (SWH) from altimeter waveform samples [2-7]. SWH is directly 
related to the second moment of the surface elevation distribution; 
there is also a preliminary altimeter measurement [5] of the surface 
skewness, related to the third moment of the surface elevation distribution. 
All of these papers use the same basic procedure: a specific probability 
distribution form (usually a simple Gaussian) is assumed for the ocean 


This work was supported by the National Aeronautics and Space Administration, 
Wallops Flight Center, under Contract NAS6-281Q. 

^ GEOS-3 had 16 waveform sample-and-hold (S&H) gates at 6.25 nanosecond 
spacing, while SEASAT-1 had 60 S&H gates at 3.125 ns spacing. 




surface, amplitude and timing biases are removed from waveform sampler 
data if necessary, a suitable waveform sample averaging time is chosen, 
and the average sampled mean return waveform is "best-fitted" by a 
theoretical template through some process of varying the template parameters 
until some typical least-squares error criterion is satisfied^. In 
several cases [3,7] differences between adjacent sampled waveform values 
are fitted to the derivative of the theoretical return waveform, but 
this is a minor modification of the basic template-fitting procedure. In 
the general template-fitting, either a simplified form is assumed for 
the different convolution terms so that the final expression will be 
simple and suited to easy computer implementation and fast computer 
running time, or a more general, more complete template is used at the 
expense of having to do several numerical convolutions for each sample 
point within the linefitting process. 

This paper presents general expressions for the radar mean return 
waveform for the case in which the ocean surface elevation distribution 
and the radar system point- target response can both be represented by a 
modified Gaussian form including skewness and kurtosis terms, and the 
antenna pattern can be assumed to be Gaussian in angle; these expressions 
give the mean return waveform at the times of interest in the region of 


‘ One exception to the template-fitting procedure is described in [8] in 
which Priester and Miller describe work intended to estimate the surface 
height probability density function without the necessity of assuming a 
specific functional form; this method is difficult to implement and 
interpret and will not be considered for purposes of this paper. 
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the risetime portion ("ramp portion") of the waveform for small off- 
nadir angles. These general expressions should be useful in either 
waveform parameter estimation from experimental data, or for waveform 
generation as part of design studies of system trade-offs for future 
radar altimeters. Comparable results can be obtained from numerical 
convolution but the expressions in this paper have the strong practical 
advantages of reducing computational time and of avoidina the detailed 
questions of sample density and round-off problems which should be part 
of any numerical convolution modeling study. 

II. EVALUATION OF MEAN RETURN WAVEFORM 

This work is based on Brown's 1977 paper [1] which reviewed the 
assumptions and limitations of the convolutional model for near-normal- 
incidence rough-surface backscattering of short-pulse radar waveforms. 

The general square-law-detected waveform W(t) is given by the convolution 

W(t) = Pps(t)*q^(t)*s^(t) (1) 

where Pp 3 (t) is the average flat-surface impulse response, q^Ct) is 
related to the surface elevation probability density of scattering 
elements ("specular points"), and s^(t) is the radar system point target 
response which includes the transmitted pulse shape. Each of the three 
terms in (1) will be separately discussed; then through a small -argument 
series expansion for Brown's Pp 5 (t), the W(t) of (1) will be expressed 
as a power series for which the first four terms are evaluated in this 
paper. 
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For cases of practical interest for satellite-carried radar altimeters, | 

the flat surface impulse response is given by Brown's equation (9) as 

Pps(t) = A exp(-6t)lQ(t^3)U(t) (2) ; 

i 

I 

in which I 

J 

6 = (4/Y)(c/h)cos 2?| (3) 

and 

3 = (4/y)(c/h)^in | 2? |_ (4) 

In (2), U(t) is a unit step function, is a modified Bessel j 

function, and 5 is the off-nadir pointing angle. In (3) and (4), c is 

the speed of light, h is the satellite altitude, and y is an antenna | 

beamwidth parameter defined as in Brown's equation (4) by a Gaussian 

approximation to the antenna gain of the form ? | 

G{e) = G^ exp[-( 2 /y)sin 20 ] ( 5 ) 

If 0 ^ is the usual antenna beamwidth, i.e., the angular full -width at 
half-power points, (3) and (4) can be written as 





and 


= luJ — (c/h)'“ sin)2e|. 0) 

sin2(0^^/2) 


The amplitude term. A, in (2) actually contains a number of different 
constants; 


A = 


G^-\-ca°{0) A 

4(4TT)=Lph3 sin-c) 


where A = radar wavelength. 


( 8 ) 


o°(0) = ocean surface backscattering cross section at normal 
incidence. 


Gq = radar antenna gain at nadir, 

and Lp = two-way propagation loss. 

In (8), the possible variation of a° with angle of incidence (or, 
equivalently, with increasing time in the return waveform) has been 
ignored; this is a reasonable assumption for sea surface scattering in 
the short-pulse (<20 nanoseconds pulsewidth ) radar altimeters considered 
here. Because radar return signals are normalized by a typical altimeter's 
AGC system, we will ignore all individual terms within the A in (8), and 
will use A as just a simple amplitude scaling term. 


B. Radar-Observed Surface Elevation Density Function 

Equation (1) includes a term q^Ct) which is the surface specular 
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point density function written in the altimeter's time domain; the 
conversion factor from surface elevation measurements in meters to two- 
way ranging time in nanoseconds is c/2 = -0.15 m/ns. The negative sign 
indicates that an increase in ranging time corresponds to a decrease in 
surface elevation^. This paper makes the usual assumption that the 
surface specular point density function is identical to the true (geometric) 
surface elevation density function, but the possibility of a difference 
between the radar-observed and the true geometric elevation densities 
must remain an open issue for future work‘d. 


This change of sign is important only for odd moments. Specifically, 
when a template-fitting procedure is used on satellite sampled waveform 
data to obtain a (time-domain) estimate of the skewness contribution 
from the surface, the desired surface elevation skewness will have the 
same magnitude but opposite sign from the time-domain-estimated skewness 
component. 

Radar measurements by Yaplee et aj_. [15] can be interpreted [16] 
as showing that the mean of the scattering distribution is displaced 
downward from the true water level by 20% of the rnis waveheight and 
recent theoretical work by Jackson [17], assuming infinitely long- 
crested waves, is in agreement with these measurements. However, 
preliminary results from the aircraft-borne Surface Contour Radar 
experiment [18] indicate that this bias may be in range of 4-8% of 
the rnis waveheight with more measurements currently underway [19]. 


In previous work the ocean surface elevation density function qg(t) 
has been assumed to be skewed Gaussian form given in the time domain 

by 


qjU) = n +~ H3(,/<,j)]expt-Ji(t/03)2], (3) 

where is a risetime, is the skewness, and is a Hemiite polynomial. 
Pierson and Mehr [9] discussed this form for radar altimeter analyses; 
this is a low order case of a general probability function by Longuet- 
Higgins [10] for a random variable that is weakly nonlinear. The first 
six Hermite polynomials for a general argument z are: 

Hl(z) = z, 

H2(z) = z— 1, 

H^(z) = z3 - 3z, 

OO) 

H^(z) = z*+ - 6 z2 + 3, 

Hg(z) = z5 - 10 z 33 + 15z, 
and ~ ~ ^ ~ 

The in (9) is related to the significant waveheight (SWH) by 

SWH = 4(c/2)a3 (11) 

in which the c/2 converts from ranging time to surface elevation, and 
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the factor of 4 is from the definition of SWH as four times the rms 
waveheight; since SWH is related to a second moment, there is no sign 
change from the time-domain-observed to the surface elevation quantity. 
Most of the radar altimeter SWH measurements to date have been based 
on the pure Gaussian resulting from setting = 0, but (9) is the general 
form used in attempts to recover a surface skewness [5]. 


Equation (9) is the result from taking the first two terms only 
in the general Gram-Charl ier series [11]. Recent work by Huang and 
Long [12], based on laboratory measurements of the surface elevation 
density function for a wind-generated wave field, suggests that (9) 
is an inadequate form for the surface elevation density, and that it is 
necessary to use the four-term series given by 


qs(t) = 


1 




o 

H3(t/„3> +^H,(t/Os) +§H,(t/Cs) 


( 12 ) 

] exp[-*2(t/ug)‘!], 


where is the kurtosis with other quantities as already defined. The 
general waveform result to be derived in the following portions of this 
paper will be based on (12), but terms will be kept separate from 
the terms so that the final results can easily be converted to a 
surface elevation density of form (9) instead of (12) if desired. 


C. Effect of Radar System Point Target Response 

The radar system point target response s^(t) is primarily the 
transmitted radar pulse shape but also includes effects of the bandwidths 
of the transmitter and receiver. When the system point target response 
is nearly Gaussian, the same general form can be used as already 
discussed for the surface elevation density. 
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1 + H3(t/a^) + 2^ H4(t/a^) + ^ Hg(t/a^) exp[-Js(t/a^)2]. 


In this case, the convolution in (1) of q 5 (t)*s^{t} cd^ itmiGdiately 
be written as 


B(t) = qg(t)*s^(t) 


a 


1 + 


H3(t/a) + ^ ^ Hg(t/0) 


exp[-*5(t/a)2]. 


in which a, X, and k (with n£ subscripts) are the composite risetime, 
skewness, and kurtosis; these composite quantities can easily be 
shown to be related to the contributing surface elevation and point 
target quantities by 


(14) 


(j2 = 0^2 + (j^2 ^ 

^ ~ ^ ^r ’ 

and K = + <^(0^/0)** . 


(15) 

(16) 
(17) 


In a practical radar altimeter, the waveform sampling gates 
will be positioned by a range tracker having its own range-noise 
(jitter) characteristics, and the probability density function of the tracker 
jitter must be included as an additional term in the convolutional 
description of the mean return waveform. For this paper the jitter will 
be assumed negligible in comparison to the other terms^. 


’’’ Usually when the jitter can not be treated as negligible, it can be 
described by a near-Gaussian form with its own 0 ., x. and k. and (15), 

J J J 

(16), and (17) will each have one additional term of the same form as the 
two terms already written. For instance, (15) will become 02 = 02 + ^ 2 + 0 . 2 . 

I> I J 
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D. Expansion and Solution for W(t) 


The waveform W(t) is given by the convolution of the PpgCt) of (2) 
with the B(t) of (14) , 


W{t) = Ppg(t)*B(t) 

CO 

(z) B(t-z)dz, 

— 


(18) 


but the presence of the Bessel function in Pp 5 (t) prevents the easy 
integration of the integral in (18). For the waveform regions of 
interest in the case of short-pulse radar altimeters in this paper, it 
is possible to expand the in a small -argument series expansion, 
given by expression 9.6.10 in [13], 


!.(z) = S (zV4)" (-^)2 
n=o 


(19) 


and this will lead to a term-by-term integration in (18) with the 
result 


where 


2 ..o 


w(t) = f exp[-d(r + d/2)] (1 )" (-^^)^c^(t) 


( 20 ) 


" ^no *''^n1 ^^^n2 


no ^ 2 , 


T r 

y'(-2)f’ 6 + xH3(z+d) 


dz 




H^(z+d) ] 


- 2^2 
e dz 


J 


( 21 ) 

( 22 ) 

(23) 
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and 


(24) 


with 




Hg(z+d) 



X = 



- d. 


(25) 


and d = fia. 


(26) 


Notice that an arbitrary time origin shift: has been written in 

(25); in template-fitting for parameter recovery from altimeter data, t^ 
is one of the parameters to be varied for the fi’al solution since in 
most altimeter designs the altitude tracker moves the waveform sampler 
set back and forth in (time) position relative to the true mean waveform. 


Carrying out the integrations indicated in (22), (23), and (24) for 
any given n will produce results in terms of G(t) and P(x), where G(t) 
is the Gaussian, G(t) = (2ir)”'^exp(-T2/2) , and P(t) is the Probability 
Integral®; two new functions of time, and E„„, are then defined by 

Cnm " D P(t) + G(t) , for m=0,l, and 2. 
nm nm nm 

These and have been worked out for n = 0,1,2, and 3, with the 
following results: 


® Often the error function erf(x) is written instead of P(x). These 
are related by erf(x) = 2P(.2'^x) - 1. P(x) can be implemented in a 
computer through the series expansions given by expression 26.2.17 
in [13 ]. 
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Dqo “ ^ ^ » 




Eqo " X(l-3d^-3dT-T^) , 

= dV4, 

EqI = (d-d3) + (|- |<,2),.dT2 . t.3^ 
= d^/12, 



5d3 d^ 
T~ - ^ 


) + (- 



15d2 

T 



+ (3d -|d3)T2 + (|-|d2)x3 


(28) 


D,^ = -3xd2 + (6 + Ad3)t, 

^10 “ ^^d + Xd3) + Xt, 

•^ll “d^ + 7 "T. 

('J- - 3d2 + ^ ) + dx + -J- , 

^ d5 ^ 

nz 2 12 ’ 

F - + dS , / 3. ^ ^3 

+ / ^ + 2 + d 3 , T ‘t 

+ r 2 ^ ^ h'^ + 2 ^: + j 


*^20 " ■'■ ®^d +Xd3) - Sxd^T + (6 + xd3)T2, 

E2Q = (“2x - 6xd2) + (6 + Xd3)x, 

D2] == (3d^ + j ) - 2d3x + •^ T^, 

E^^ = (-2d - 2d3) + (-1 +|‘*)r. 


(29) 


(30) 
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(30 cont'd. ) 


^22 ~ + ^ ) “ d^T + Y2 » 

^22 “ ~ ’*' t^2^^ ~ ^ ’ 

D^q = (-6A - 9\d^) + (18 + 18\d + 3Xd^)T - 9Xd^r^ + (6 + Xd^)t^, 

E 3 Q = (12 + 18xd + 2xd3) - 9xd"-t + (6 + Xd3)^^ 

83 ^ = (- 6 d - 3d3) + (9d2 + _ 3 j 3^2 + d"^ 3 ^ 

En, = (4 + 9d^ + 7 ) - 3d3t + T t 2 , 

^ ^ (31) 

D 32 = (-10d3 - |dS) + i^d^ + ^ 13 , 

E32 = (-^ + ^d^ + 10d‘* )i + ( ^ + 
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III. SUMMARY AND EXAMPLES 


In brief review, the waveform is given by (20), with the number of 
terms required being a function of time, pointing angle, and antenna 
beamwidth'^. For each term C^(t), (21) gives C^(t) in terms of c^q, 

C^-j , and C ^2 with these three quantities in turn given by (27) which 
uses (28), (29), and (30), or (31) for n=0,l,2, or 3, respectively. This 
procedure may appear cumbersome and it does assume the use of a computer, 
but these results are very much easier to use and produce answers in far 
less computer time than the numerical convolutions that would otherwise 
be needed®. 


For instance, a single term is adequate fov the GEOS-3 altimeter 
waveforms at any time within the span covered by the individual sampling 
gates for pointing angles within 1° of nadir, while SEASAT-1 with its 
smaller beamwidth and larger sampling gate span in time requires the 
first three terms in (20) for any position within the sampling gate span 
and for pointing angles within 1° of nadir. 

® One of the uses of the waveform expansion (20) is as the template to 
fit to radar altimeter data. Typical non-linear iterative least-squares 
linefitting routines require the derivatives of the waveform expression 
with respect to the individual parameters as evaluated at each input 
sample point. Such derivatives may be estimated by numerical differences, 
varying each parameter concerned. Alternatively, the derivatives may be 
obtained by differentiating (20). While the resulting expressions are 
probably not of sufficient general interest to justify their appearing 
in this paper, they have been worked out and checked and are being used 
in linefitting routines by the author; these derivative expressions will be 
mailed upon receipt of a request. 
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Figures 1-6 provide examples of mean return waveforms calculated 

from (20) for a radar altimeter having an antenna beamwidth 0 ., = 1.6 

w 

degrees at a height 8 x 10^ meters above the ocean and having a pure 
Gaussian radar system point target response of 3.125 nanoseconds full- 
width at ’2 height so that = 1.327 ns, and = 0 in (14). These 

numbers are nominal SEASAT-1 values £14], and these Figures predict 
SEASAT-1 measured mean return waveforms for the limit of infinite averaging 
time^. 


The actual SEASAT-1 mean return waveforms will differ from this 
paper's Figures because of the distinctly non-Gaussian sidelobe structure 
in ^he transmitted pulse shape. MacArthur, on pg. 4-50 of [14], describes 
an attempt to treat this sidelobe structure by increasing the effective 
width of the Gaussian used to represent the transmitted pulse; however, 
a discussion of the full consequences of the actual SEASAT-1 transmitted 
pulse shape is beyond the scope of this paper. 
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Figure 2. Mean Return Waveforms in Idealized SEASAT-1 Radar Altimeter 
at 0.4® Pointing Angle, Showing Effects of Number of 
Terms in Waveform Expansion. 
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Figure 3. Mean Return Waveform in Idealized SEASAT-1 Radar Altimeter 
at 6.0 m SWH, Showing Effects of Different Off-Nadir 
Pointing Angle 
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Figure 4. Mean Return Waveforms in Idealized SEASAT-1 Radar Altimeter 
at 0.0° Pointing Angle, Showing Effects of Different Ocean 
Significant Waveheight 
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